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ABSTRACT 

We present the stellar velocity dispersion measurements for 5 Luminous Compact 
Galaxies (LCGs) at z=0.5-0.7. These galaxies are vigorously forming stars with average 
SFR ~ 40 Mq/ji. We find that their velocity dispersions range from ~ 137 km/s to 
260 km/s, while their stellar masses range between 4 x 10^ and 10^^ Mq. If these LCGs 
evolve passively after this major burst of star formation, their masses and velocity 
dispersions, as well as their evolved colours and luminosities are most consistent with 
the values characteristic of early-type spiral galaxies today. 



1 INTRODUCTION 



Luminous Compact Galaxies (|Hammer et al.l l200l[ ) are 
starbursts galaxies at intermediate redshifts mostly de- 
tected in both UV and IR wavelengths, characterised by 
having small effective radii {Re < 5 kpc), high lumi- 
nosities (Mab ( B) ^ —20) and strong emission lines. 
iHammer et al.l (|200ll ) observed a representative sample of 
LCGs selected from the CFRS fields using the interme- 
diate resolution {R > 600) spectrograph FORSl and 
F0RS2 on the VLT/Kuyen telescope. The spectra revealed 
some strong absorption lines (Ca II K and H, G Band, 
Fe I, and Balmer lines) as well as narrow and intense 
emission lines ([OII]A3727A, [OIII]AA4858, 5007AA, Balmer 
lines ) . The spectro-phot o metric analy sis of these galax- 
ies (jHammer et al.l l200ll . iGruell I2OO2I ) showed that they 
are likely composed of three different stellar populations. 
The youngest population presents strong emission lines 
([OIIJA3727A Balmer lines), indicating present day active 
star formation and sub-solar metallicity. A second stellar 
population was formed within the last few hundred million 
years. It is characterised by a solar metallicity and the pres- 
ence of Balmer lines in absorption. The third stellar pop- 
ulation, older than 5 Gyr, exhibits solar metal abundances 
and strong metallic ab sorption line s (Calcium, Iron lines, G 
band, Titanium, etc.) (|Gruelll2002l ). 

Almost all LCGs analysed by iHammer et all (|200ll ') 
have large extinction coefficients {Av ~ 1-5), yielding 
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average extinction-corrected star formation rates SFR ~ 
4OM0 yr~^ (which is ~10 times higher t han t h ose estimate d 
from the UV fluxes) ((Hammer et ali boOll . ICruell I2OO2I ). 
Some LCGs show morphological irreg ularities and/or close 
companions as rev ealed in HST images l|Hammer et al.|[200ll . 
IZhenget al]|2004h . These observations suggest that LCGs 
may be undergoing violent star formation e vents similar to 
those occurring in close interacting systems. IHammer et al.l 
(|200ll ) concluded that these LCGs may be the progenitors 
of the bulges of spirals galaxies forming inside-out. 

A new parameter to shed hght on LCGs local counter- 
parts is the galaxy's kinematics. Measurements of the Hp 
emission line velocity widths for Hammer et al. LCG sample 
showed low velocity widths of ~ 70 km~^, suggesting they 
may be low mass objects (~ lO^M©). Ha velocity widths 
for a handfiil LCG s measured with ISAAC at the VLT 
(|Tresse et al.ll2002l ). showed a "double horn" profile char- 
acteristic of rotation. Integral spectroscopy was also used to 
measure the velocity field for some LCGs by iPuech et al.l 
(.2006 ) using GIRAFFE at the VLT. However, due to the 
small apparent size of these objects, the velocity map ex- 
tends only over very few spaxels. 

The most reliable measurement of galaxy kinematics is 
the stellar velocity dispersion from absorption lines. In this 
paper, we present the first velocity dispersion measurements 
from absorption lines for 5 LCGs at intermediate redshift. 
We constrain our study to the spectral range that includes 
the Balmer lines {HW, H9, H8, He, H5, H-/), the calcium 
doublet Call K,H and the G Band. Section [2] presents the 
data set for our galaxy sample. Section[3]describes the meth- 
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ods used to measure the velocity dispersion, the photometry 
and the stellar masses of these galaxies. Section |4] contains 
the results and the discussion. We assume the following cos- 
mology in this paper: Ho = 70 km sT^ Mpc~^, Q,\ = 0.7 
and Q,m ~ 0.3. 



2 DATA 

2.1 Sample selection 

The galaxy sample was selected from three Canada 
France Redshif t Surv ey (CFRS) fields: CFRS 0000+00 
llLe Fevre et all 1 19951 ). CFRS 0300+00 l|Hammer et all 
1 19951) and CFRS 2230+00 jLillv et al.lll995l ). Intermediate-z 
LCGs were selected using criteria defined by [Ham mer et al. 
lj200ll ): size (ri/a ^ 5h^o kpc), luminosity (Mab(B) ^ 
-20.0), redshift (0.5 2 1) and the pres- 

ence of a major star formation episode characterised 
by an [OH]A3727A emission with an equivalent width 
EW([OII]X3727k) > I5A. This resuhs in a sample of 32 
LCGs, or 29% of the most luminous (Mas(B) ^ -20.0) 
galaxies found in the three CFRS fields (Gruel 2002). We 
observed 22 of these galaxies with the FORS/R600 and 1600 
spectrograph at the European Southern O bservatory 8m 
VLT/Kueyen at a resolution of FWHM=5 A ([Hammer et al.1 
I2OOII . lGrueill2002l '). The typical exposure times were 12000 
seconds per object. 

We reduced the spectroscopic data by using the Mul- 
tired package within IRAF[3 Measurements and analysis 
(flux, equivalent width, star absorption co rrection, met allic- 
ity) were done with software described in lCruell ()2002l ). 

In this paper, we analyse a subsample of give LCGs with 
the higest S/N and the strongest absorption lines, for which 
measurements of stellar velocity dispersion are feasible. 



3 METHODS 

3.1 Velocity dispersion measurements from the 
absorption lines 

Velocity dispersion was measured with the program 
Movel, included in the REDUCEME softwares package 
(jCardiel et al.l 1 19981 ). The instrumental resolution of our 
LCG sample is FWHM~ 5 A (cri„s ~ 100 km/s). The 
lowest velocity dispersion measurable with this instrumen- 
tal resolution is cr ~ 70 km /s for spectra with S/N > 10 
iMatkovic fc Guzman|[2005l ') . We thus selected galaxies from 
our LCGs sub-sample with S/N > 10 per resolution ele- 
ment. The S/N criteria reduced the LCGs sample to only 5 
objects. The final sample is given in table [2l 

Since the original aim of the observations was to anal- 
yse emission lines, we did not observe stars that we could 
use as templates. Therefore, we created our own a series of 
stellar templates with the same instrumental resolution as 
the LCGs. 

We used stellar templates from MILES spectral stellar 

^ IRAF is distributed by National Optical Astronomical Obser- 
vatory which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Example of one stellar spectrum (22.0637) (plain line) 
with the best broadened stellar template fit (dashed line) obtained 
with Movel. The different absorption lines used are labelled. 



library (jSanchez-Blazauez et al.ll2006l . iFalcon-Barroso et al.l 
201lj, which includes stars observed between July 2001 and 
December 2002 in La Palma (Spain) at the Isaac New- 
ton telescope. These templates are the basis of our galaxy 
template. Since their resolution is higher (~ 2. 5 A FWHM) 
than our sample of galaxies (R=600), we convolved them 
with a gaussian kernel. The characteristics of the gaussians 
used for the convolution was determined from the differ- 
ence between the instrumental resolution of the galaxies 
and the stellar database. The instrumental/spectral reso- 
lution of the individual LCGs was determined from the sky 
spectrum obtained with the same slit as the spectrum of the 
galaxy. The intrinsic dispersion of the instrument is given by 
f/nst = ^^2^35°''^ X i X c. The valuc of the Gaussian func- 
tion used to broaden the stellar template is the quadratic 
difference of the instrumental resolution of our observations 
and that of the stellar library: Act = <^jnst ~ '^Itars (See 
Fig.dJ. 

Using a stellar database to build stellar templates by 
convolving the instrumental resolution to one's own observa- 
tions has already been succes sf ully d er nonstrated f or early - 
type galax i es by iTreu et al.l l| 19991 ), iTreu et al.l (|2001al ). 
iTreu et al.1 l|2001bl ). The studies proved that at any given 
S/N and resolution, there is a lower limit to the velocity dis- 
persion that can be measured. Typically, this limit stands at 
half the inst rumental resolution of a given ga laxy spectrum 
at low S/N l|Bender. Burstein. fc Faberlll992h . All measure- 
ments with a S/N < 10 per reso lutio n element are deemed 
not reliable. iTreu et al.1 (|2001ah and iMatkovic fc GuzmanI 
( |200f^ quantified the errors for different S/N and range of 
velocity and showed that the systematic error for a galaxy 
with a S/N> 10 and a velocity dispersion of 150 km/s is 
~ 15%. 

Velocity dispersions for the galaxies were measured us- 
ing the MOVEL and OPTEMA algorithms described by 
Gonzalez (1993). The MOVEL algorithm is an iterative pro- 
cedure based in the Fourier Quotient method (Sargent et al. 
1977) in which a galaxy model is processed in parallel to 
the galaxy spectrum. The main improvement of the proce- 
dure is introduced through the OPTEMA algorithm, which 
is able to overcome the typical template mismatch problem 
by constructing for each galaxy an optimal template as a lin- 
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ear combination of stellar spectra of different spectral types 
and luminosity classes. 

For this study, error in velocity dispersion measurement 
was statistically determined by bootstrapping. One hundred 
pseudo-spectra were randomly simulated in the range given 
by the galaxy error spectrum using a gaussian noise model. 
The error of the velocity dispersion measurement was ob- 
tained from the statistical distribution of the 100 pseudo- 
spectra measurements. 

Different effects limit the precision or the ability to mea- 
sure galaxy velocity dispersion. Stellar template mismatch is 
one such limiting factor. This effect happens when the star 
template used to fit the galaxy spectrum are not represen- 
tative of the real stellar population of the galaxy. Creating 
a galaxy template from a stellar library can be impossible, 
the stellar library does not contains every stellar type and 
one of the main stellar population which composed the ob- 
served galaxy can be missing. As this template is a first guess 
to start the velocity dispersion measurement, no numerical 
criteria to avoid template mismatch is used. It is possible 
to help the software to make this first guess by removing 
some extra-features (see below) or limiting the number of 
stars in the library to the one presenting the most impor- 
tant features for velocity dispersion measurement. We re- 
moved every galaxy spectrum from our sample when an im- 
portant template mismatch at the initial step was visually 
detected even if their S/N was over 10. To test the impact of 
a small initial template mismatch (no visual detection and 
X'^ greater than the template automatically created by the 
software) we forced the software to use specific templates, 
thus creating an artificial template mismatch. Measurement 
of the velocity dispersion with this mismatch showed an er- 
ror of ~ 5% in the dispersion measurement. 

In our case, another effect is introduced by the convolu- 
tion of the instrumental resolution of the stellar templates to 
that of the LCG's. The effect of this modification was inten- 
sively tested with galaxies of well known velocity dispersion. 
An elliptical galaxy from the Coma cluster was changed and 
the velocity dispersion measured with the templates con- 
volved to a slightly different one. The procedure was re- 
peated with a pseudo-instrumental resolution modified up 
to ±3A in O.O5A steps. An error of O.lA in the difference of 
resolution introduces an error of 5% in the velocity disper- 
sion measurement. The error in velocity dispersion increases 
more strongly with higher error in resolution. At ±0.1 A the 
error is ~ 5%, while at ±0.2A and ±0.5A the errors are 
respectively ~ 10% and ~ 22% in the dispersion measure- 
ment. The galaxies instrumental resolution were measured 
using different lines from the sky spectrum extracted at the 
same t i me th an the LCG spectra with the same polynomial 
l|Gruell 12002! ) . The resolution measured from the different 
lines was stable at ±0.1 A along the spectrum. 

Small emission lines tend to fill the Balmer lines (our 
major features) changing the global shape of the absorption 
features To quantify the infiuence of this contamination, the 
emission lines in our sample were masked and velocity dis- 
persion was measured without them. We also observed that 
if not removed entirely, an artifact can appear at the contin- 
uum substraction step for the strongest emission lines such 
as [OIIJA3727A or H/3 and bad sky substraction features. As 
a result these lines were erased manually. Both effects were 
quantified (measured with and without the lines) as an er- 



ror of 2% for the strong emission line [0//]A3727A and 5% 
for the emission inside the Balmer lines. For galaxy 03.0645, 
where the absorption lines were weaker and contaminated 
by inside emission, inside emission lines were also removed 
and led to an increased error in velocity dispersion. 

To emphasise the strongest or the most useful ab- 
sorption lines, the analysis has been restrained to the 
wavelength range AA[3570, 4400]AA. The best absorptions 
lines present in the spectra were used (Balmer lines 
(HIO, m, H8, He, HS, Hj), Call K and H and the G Band). 
These restriction have no incidence in the measurement and 
improve continuum substraction. 

As a summary, the error analysis showed that the accu- 
racy of velocity dispersion measurement is principally lim- 
ited by the noise in each galaxy spectrum. Our LCGs sam- 
ple was thus restrained to galaxies with S/N > 10. The 
second limiting effect is the possible mismatch in resolu- 
tion between the stellar templates and the galaxies. This 
effect was minimised by adjusting the stellar resolution to 
the galaxy spectra. Still an error of 10% was found for a 
typical error of O.lA. The average template mismatch was 
evaluated at around 5%. Finally error due to the presence of 
strong emissions lines, an error purely numerical, was found 
to have minimal impact with value of less than 2%. 



3.2 Photometry 

The B, V, I and K photome try for our sample were obtained 
from the CFRS catalogues l|Lillv et al.ll 19951 . iLe Fevre et~al] 
Il995l . iHammer et al.1 119951 ). The absolute AB magnitudes 
and rest-frame colours were derive d from the best fit to 
the photometries data of models from Bruzu al A. fc Chariot 
(1993) and Br uzual fc C hariot (2003). The errors were de- 
rived using monte-carlo simulations and estimated to be 
~ 0.14 and ~ 0.2 mag for the magnitudes and colours re- 
spectively. Note that all the values were transformed to the 
concordance cosmology used in this paper. 



3.3 Stellar mass measurements 

Stellar mass were calcul a ted u sing t he code descr i bed in 
ICristobal-Hornillos et akl ll2005h and iHempel et al] (|201ll ) 
based on Guzman et al. (|2003l ) idea. The ages and masses 
were determined by fitting the observed photometry (B, V, 
I and K band) to the modelled fiux obtained from the con- 
volution of a redshifted library of synthetic galaxy spectra 
with the filter transmission functions. 

A two component synthetic model consisting on a young 
instantaneous burst and an exponentially declining SFR was 
considered. Models were built using B C03 stellar population 
predictions (|Bruzual fc Chariot! bOOlj ) considering different 
ages for both components as indicated in Table[T] The mod- 
els depend on several parameters: IMF, metallicity, Av, ex- 
tinction law, mass, SFR, and age. We considered the same 
IMF, metallicity, Av and extinction law for both old and 
young component. Therefore a total of 9 parameters (given 
in TablelTJ define the two-component population model. To 
handle this degenerate pro blem, we first made simulations 
(|Crist6bal-Hornillos! l2005bf ) to understand the infiuencing 
parameters on the stellar mass determination. We consid- 
ered a modelled LCBG galaxy with the mean parameters of 
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iGuzman et all dioosl ) and placed it at z = 0.5, z = 0.8, and 
z = 1.1, adding realistic photometric uncertainties (0.03, 
0.06, and 0.1 mag for each redshift). For the simulations we 
considered the U, B, V, I, K bands. We saw that extinction 
and metallicity were the most influencing parameters for fi- 
nal stellar mass. Increasing the amount of extinction tended 
to decrease the inferred stellar mass. Due to the extinction- 
age degeneracy, a younger and less massive underlying com- 
ponent increased the extincted flux in the rest-frame UV 
bands. Similarly, a lower metallicity also increased the stel- 
lar masses. Other parameters such as SFR, extinction law 
and IMF had little effect on the estimates of masses. In our 
study, the best models (x^ ^ Xmin + 1) deviated less than a 
factor 2 in the inferred stellar masses. 

We fixed certain parameters using measurements from 
the real spectra. The extinction Av was measured from the 
Balmer decrement between II/3, and Ha emissions lines 
and infrared flux when available. In the initial approach we 
used the Large Magellanic Cloud extinction law for all the 
galaxies. For one of them, 03.1540, the models could not fit 
the SED correctly. The effect was reduced by using Calzetti 
law for the extinction law but still, the fit did not reproduce 
the expected SED. The mass determined for this galaxy was 
at the upper limit and probably overestimated. This LCG 
was also detected by ISO and considered peculiar in both 
term of its extinction and star formation. 

For the fitting process, the ranges considered for the 
different parameters are shown in Table[T] The is cal- 
culated from the difference between observed and modelled 
photometry as in equation[T] 



i 



rpobs 
j 



(1) 



where the weights aj in each filter, are given by 
1.086/a"mj where is the error in magnitudes. At each 
iteration, once the model parameters are fixed, the fitting 
process determines which combination of stellar masses and 
ages for the components produces the best fit to the ob- 
served photometry. A simultaneous fit to the two-component 
model is performed. We attempted an alternative approach, 
consisting of first fitting the young burst component to the 
blue bands photometry, and then fitting the residuals to 
the underlying component. Such approach seemed reason- 
able since the bluest /reddest bands were expected to be 
dominated by the young/old component. However, the re- 
sults were found to strongly depend on which of the two 
components was fitted first, i.e. if we first fitted the young 
component, then too much NIR fiux was assigned to the 
component, yielding a low total mass. 

The upper age was defined as the age a galaxy formed 
a,t z — 4 would have at the z observed using our adopted 
cosmology. We took 10^ yr as the lower age limit for re- 
cent burst of star formation, since the models do not in- 
clude nebular component. The lifetime of the stars produc- 
ing most of the ionising photons are typically 3 — 5 x lO" yr; 
after 10^ yr t he production rate of ionis ing photons drops 
by over 99% JCharlot fc Longhettill200ll '). The simulations 
(jCristobal-Hornillod l20o'5bF showed that in a two compo- 
nent fitting model, when the burst is allowed to be younger, 
the underlying component turns out to be older and more 



Table 1. Ranges of the parameters used to model fitting. 



Parameter 



Range 



Age of recent burst 10^ — 5 X 10** yr 

Age of underlying population ... 2 X 10* — 2 X lO^'' yr 

IMF (Salpeter 19551 

Mass of stars: lower, upper limits 0.1, 100 Mq 

Av Fixed 

Extinction Law Fixed 

Metallicity 0.4 Z©, Z© 

SFR (underlying) r = 1.0Gyr 

SFR (burst) Instantaneous 




10000 15000 

A(A) 



Figure 2. Fit of one LCG (22.0637) photometry with two dif- 
ferent galaxy populations. Plain line, the model is fitting the ob- 
served photometric points and the dashed line is the model which 
has passively evolved until z=0. The evolved photometric point 
are measured by integrating this spectrum. 



massive. In that sense the stellar masses estimated here are 
a lower limit. 

The masses measured for the LCGs sample using our 
own code indicate that these galaxies have a stellar mass 
4.57 X IO^Mq ^ M* < 1.38 x IO^Mq. 

To estimate the evolution of these galaxies to z — 0, 
the code fits the galaxy at the observed redshift and we 
considered the fitted models as they continued to evolve until 
2: ~ 0. We assumed that all gas in the fitted galaxy had been 
used to produce stars. 

In the models, we are supposing a passive evolution for 
our LCG. They are experiencing their last starburst and will 
not have any mergers. This last hypothe sis is sustained by 
IConselice. Blackburne. fc PapovichI |2005| who showed that 
the major merger rating decrease since z ^ 1. Under these 
conditions, velocity dispersion and masses are independent 
of the evolution of the galaxies and the colour and luminosity 
will be determined by single evolution of the best fit of the 
observed SED. The evolved point used later in this paper 
was calculated from the models after evolution (see Fig. [2]). 
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Figure 3. In these figures, the points are the LCGs analyse in 
this paper, with a triangle for the ISO galaxy not very well fit- 
ted by D. Cristobal models. The grey points are the position on 
the LCGs at the time of the observations and in white the same 
objects after passive evolution to 2; = 0. To find the nature of 
the LCGs we divided the SDSS data, the reference sample used 
here, in three different categories: Ellipticals, Spirals and Irregu- 
lars g alaxies. The crosses repre sent a local sample of HII galaxies 
from l|Telles fc Terlev ich"l997^. The upper row is the M* versus 
a diagram, the middle one the Faber-Jackson relation L versus a 
and the lower the colour versus a relation. 

4 RESULTS AND DISCUSSION 

The velocity dispersion measurements and the stellar mass 
estimates are listed in the table [2l The velocity disper- 
sions from the absorptions lines of the LCGs range between 
~ 137 km/s and 260 km/s with a median of ~ 180 km/s. 
The two galaxies with highest velocity disper sion (03.1349 
and 03.1540) are ISO galaxies. [Hamme 3 1 19991 ) showed that 
galaxies detected by ISO tend to be large and massive, and 
found mostly in interacting systems yielding strong star for- 
mation rates (> lOOM© yr~^). 

In FigureO we compare stellar masses, Mb and iJJ — V) 
colour of LCGs with different type of nearby galaxies. 

The data used for comparison include several different 
kinds of galaxies taken from the SDSS stellar mass cata- 
logue (Kauffmann et al.ll2003l ). We remark that this cata- 
logue does have some degree of mismatch within a 1" ra- 
dius with the DR4 catalogue we used for the magnitude 
and colours. We removed all objects whose cross correla- 
tion yielded different plate identification number, different 
fibre identification number, different modified Julian Date, 
discrepancies in the redshift or a large error in the mea- 
surement of the mass or the velocity dispersion. Our final 
selection include ~94% of the initial sample. 

Figure [3] summarise the results of our analysis. The 
data points are the 5 LCGs studied in this paper. In grey, 
the observed value, in white the model predictions after 
passive evolution. The triangle corresponds to the object 
whose Stellar Energy Distribution (SED) determined by our 
stellar mass code did not fit its photometry values. There- 
fore its estimated evolution is very uncertain. The crosses 
in the colour-cr diagram are HII galaxies (see figure cap- 
tion). The grey scale are the density plots for different SDSS 



galaxies. The separation in different galaxies types shown in 
the three columns was done by colour selection, following 
' Pukugita. Shimasaku, fc Ichikawal (liggsl ). 

The upper row is the M* versus a diagram, the middle 
one the L versus a diagram, and the lower the colour versus 
a diagram. About ~5 Gyrs ago, the 5 LCGs have stellar 
masses and velocity dispersions consistent with those char- 
acteristic of today's population of early-type spheroids and 
spiral galaxies. Their luminosities and colours, however, are 
similar to most luminous, bluest local HII galaxies. Assum- 
ing the stellar masses and velocity dispersions of these LCGs 
remain approximately constant, and that they are undergo- 
ing their last burst of star formation, simple evolutionary 
synthesis models predict that these objects will evolve pas- 
sively to best resemble the typical luminosities and colours of 
early-type spiral galaxies. Our velocity dispersion and stellar 
mass measurements, combined with this simple evolutionary 
predictions, are consistent with th e proposed link betw een 
LCGs and massive spiral galaxies (|Hammer et al.ir200ll ). 

We note that our results are not necessarily inconsis- 
tent with previous works suggesting that other class of in- 
termediate redshift star-forming galaxies (the so-called Lu- 
minous Compact Blue Galaxies) may evo lve into today's 
population of low mass spheroidal galaxi es (|Koo et al.ll2005l . 
iGuzman et "aLll2003l . iNoeske et al.ll2006l ). The LCGs in our 
sample have a luminosity similar to the 15% brigh t est ob - 
ject s in the LCBGs sam ples studied by IKoo et al.l l|2005l ), 
and lNoeske et all (|2006l ). 
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S/N 


z 


M{B)ab 


A(V) 


[O/H] 


^[O/H] 


<^abs 




log{M^) 








(mag) 


(mag) 






(km/s) 


(km/s) 


(Mq) 


03.0645 


12 


0.527 


-20.45 


1.53 


8.57 


±0.09 


133.3 


±91.2 


10.68 


03.1349 


13 


0.616 


-21.15 


1.08 


8.98 


±0.13 


196.5 


±15.9 


11.14 


03.1540 


12 


0.689 


-21.27 


3.52 






259.7 


±37.3 


10.89 


22.0429 


12 


0.624 


-20.02 
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